Hardware-based Demonstration of Time-Delay Interferometry 
and TDI-Ranging with Spacecraft Motion Effects 
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This hardware-based experimental simulation of the Laser Interferometer Space Antenna (LISA) 
implements a real-time electronic duplication of the time-changing inter-spacecraft (SC) laser phase 
delays while measuring heterodyned laser fields with /icycle phasemeters. The beatnote measure- 
ments are used to verify the capabilities of theorized post-processing time-delay interferometry 
(TDI) combinations in the proper time-scaled and time-shifted linear combinations. The experi- 
ments meet the 18 pm/ 'v/Hz LISA measurement sensitivity goal after demonstrating the cancellation 
of 100 Hz/v^Hz laser phase noise by TDI-ranging the time- varying arm-length to an accuracy better 
than 2.0 m using a frequency modulated ranging tone. 

PACS numbers: 04.80.Nn, 95.55.Ym, 07.87.+V, 07.60.Ly, 42.87.Bg, 07.05. Hd, 42.55.Xi, 42.60.Fc, 42.60.Mi 
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INTRODUCTION 

Future space-based gravitational wave (GW) interfer- 
ometers [l[, such as the Laser Interferometer Space An- 
tenna (LISA) 0,13, will measure gravitational radiation 
from compact-star binaries and binary black hole merg- 
ers in the 0.1 mHz to 0.1 Hz frequency range, providing a 
new window through which to observe these astrophys- 
ical systems LISA exploits a modified Michelson 

GW detection technique by taking one-way laser phase 
measurements between laser benches on three individual 
spacecraft (SC) (Fig.Q] 5]) to measure and extract the 
GW spacetime strain. The SC, defining the vertices of a 
triangular constellation, follow independent heliocentric 
orbits resulting in unequal, time-changing interferome- 
ter arm-lengths. Thus, the GW measurement sensitivity 
depends heavily on the ability to combine these individ- 
ual SC data-sets to account for the laser phase noise, 
clock noise, and spacecraft motion. These time-scaled 
and time-shifted linear combinations, referred to as time 
delay interferometry (TDI) combinations 0], complete 
the laser transfer chain, cancel the laser phase noise, and 
extract the GW strain information. LISA Simulator 0], 
Synthetic LISA [3, and LISA Tools 0] have produced 
numeric simulations of these data-sets for mock LISA 
data challenges (MLDCs) [10(. Hardware based labo- 
ratory experiments have also verified Sagnac-type TDI- 
combinations with clock noise corrections using a 1 meter 
test-bed. 
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Taking the next step in validating the effectiveness 
of the TDI combinations, the University of Florida has 
constructed a hardware-based LISA simulator which uti- 
lizes a real-time digital electronic replication of the multi- 
second laser phase delays between individual laser bench- 
tops while simultaneously taking low-frequency phaseme- 
ter (PM) measurements of LISA-like photodetector (PD) 
beatnotes Previous experiments have generated 
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data-sets and tested the capabilities of the TDI-Xi com- 
binations [13 which cancel the laser phase noise in a 
static interferometer. Advancements to the simulator 
have provided the ability to model time-changing delays 
and incorporate the SC-motion induced laser phase cou- 
pling into the measurements. The following experiments 
use the TDI-AT2 combinations to cancel ~ lOOHz/yTIz 
laser frequency noise in a LISA-like interferometer re- 
sulting in greater than 10 orders of magnitude noise sup- 
pression below 1 mHz and meeting the 18 pm/\/Hz LISA 
measurement sensitivity goal in a majority of LISA-like 
experiments. The analysis also shows how the time- 
varying inter- SC round-trip arm-lengths are ranged to 
an accuracy of < 2.0 meters utilizing a TDI-ranging 
reference tone. 



MODELING LISA 

The sensitivity of the LISA detector is determined 
through a combination of requirements [3 which are de- 
fined to optimally observe scientifically interesting astro- 
physical sources while staying within the bounds of cost 
and feasibility. Each element in the LISA measurement 
chain has a pre-defined requirement (Table[TJ in order to 
meet the overall measurement sensitivity. Based on these 
specifications, the dominant sensitivity-limiting terms in 
the LISA design are the disturbance reduction system's 
(DRS) acceleration noise at low frequencies, / < 3 mHz, 
and the interferometry measurement system's (IMS) dis- 
placement sensitivity at high frequencies, / > 3 mHz. 
The DRS is implemented to limit non-gravitational ac- 
celeration noise on the six, gravitationally sensitive, proof 
masses. This will be verified with a pre-LISA test mis- 
sion, the LISA-Pathfinder [13. The IMS is responsible 
for measuring the one-way differential length between 
these proof masses. 
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TABLE I. LISA Characteristics & Requirements Q 
Characteristic Specification 



Laser Stab. 
PM Precision 
IMS Strain Sens 




DRS Accel. Noise + (g^j + gj 

Ranging Accuracy SL = 1 meter, <5r = 3.3 ns 

Arm-Length L = 5.0 ± 0.1 Gm 

Light-Travel Delay r = 16.67 ± 0.33 s 

Relative Velocity v — ±20 m/s, /3 = ±66 ns/s 

= O.lmHz, /m = 2.8 mHz, f H = 8mHz 
a Note that LISA, as a combined NASA/ESA funded mission, no 
longer exists and has since been replaced by NGO/eLISA in 
Europe while NASA develops new space-based interferometer 
concepts under the acronym SGO. However, most space-based 
GW interferometers will be generally LISA-like with similar 
measurement complications. I . I j jj The requirements 
specified and experiments performed in this description focus 
on an approximate 'worst-case' LISA-like scenario. 
b The IMS strain sensitivity refers to a single link requirement 
including shot-noise, path-length noise, residual laser phase 
noise, phasemeter noise, and many other technical noise 
sources. The experimentally relevant terms are used to 
calculate the TDI-X2 displacement equivalent sensitivity goal. 




FIG. 1. Model of the LISA-IMS: The pre-stabilized lasers, 
proof-masses, optical benches, and inter-SC light field trans- 
fers of the LISA mission are illustrated. The spacecraft, Sd, 
the associated light-travel time delays between the SC, T q / q i, 
and the photodetector observables, s sr , b sr , and, f sr , are la- 
beled. 



The LISA intcrfcrometry measurement scheme consists 
of six nearly-identical laser benchtops, two on each of the 
three SC. Each benchtop includes a pre-stabilized laser 
source, an optical bench, a proof-mass, a DRS, a fiber 
coupler to transfer the laser field between adjacent bench- 
tops, and a telescope to transmit the laser field to the 
adjacent SC (Fig.[J) @. Each optical bench uses /ucycle 
PMs to measure the differential laser phase of hetero- 
dyned laser fields on three primary PDs. The measured 
signals are encoded with the local-SC to far-SC distance, 
s sr , the local-SC to local-proof-mass distance, b sr , and 
the differential laser phase induced by the fiber back-link 
between adjacent benchtops, f sr . The TDI combinations 
of these observables are derived to complete the phase 
transfer chain and cancel the dominant laser phase noise. 

To clarify the TDI analysis we consider two simplifica- 
tions to the system. First, we assume that the counter- 
propagating fiber back- link terms, f sr , can be measured 
and accounted for and, thus, we may treat the SC as hav- 
ing only one laser source [19(. Also, assuming the DRS 
system works well enough to shield the proof-masses from 
non-gravitational acceleration noise sources and that the 
local SC to local proof-mass distance can be accurately 
measured and accounted for, it is then justifiable to inter- 
pret the SC themselves as the interferometric GW proof 
masses and the b sr terms can be neglected in the analysis 



as well [18] . This leaves the interesting TDI terms, 



S sr 



•(*) 



'Js;q 



Ma q (t-T q (0)))+g q (t), (1) 

■ g q , (2) 



or effectively, a comparative one-way measure of the local 
SC to far SC distance with a first-order velocity correc- 
tion. In this notation <p s is the phase of the 'sending' 
laser (from the adjacent SC), <fi r is the phase of the 're- 
ceiving' laser (on the local SC), and g q is the GW mod- 
ulated laser phase on the arm opposite SC q . The light- 
travel time-delays between the space-craft can written as 
r q (t) = r q (0) ± (3 q t where r g (0) is the initial light-travel 
time-delay along the arm opposite SC q , from SC S to 
SC r . fi q — [1— a q ] — v q /cwhei:ev q is the differential posi- 
tion, or velocity, between SC r and SC S . The colon nota- 
tion is used to transfer laser fields between moving frames 
by taking the time transformation, t — > a q (t — r q (0)), as 
shown in Eq.[21 of which can be applied successively as 
in Eq.2] [2(|. Note that the counter-propagating inter- 
SC light-travel time-delays along a single arm are not 
equal (t 9 (0) 7^ T q >(0)) due to the orbital rotation of the 
constellation, although in this analysis, the first deriva- 
tive is: /3 q = /3 q > [2l[. The prime notation refers to the 
different outgoing (un-primed) and returning (primed) 
light-travel-time laser phase delays. It is assumed that 
\v q \ < 300 m/s, or, |/3| < 10" 6 ; thus, 2 nd order special 
relativity corrections of order /3 2 are negligible 2^, 2^] ■ 
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TDI Theory 

Choosing a master SC, SC\, as the interferometer ver- 
tex (beam-splitter), we can cancel the delayed SC2/3 
laser phase terms from the local, s s i, signals and con- 
struct two differential, round-trip single-arm measure- 
ments by forming, 



A s =s sl + s ls - q/ 

A s =<j)l - <j) s -. q > + (<j) s - 4>l;q);q> 

A s =(j)l - 4> s;q ' + 4> s . ql - 4>i m > 

A s =4>i - 4>i m > 



(3) 



(4) 



or, explicitly as a function of time: 

A s {t) =s al (t) + s la {a g ,(t- v(0))) (5) 
A s (t) =<^{t) + Ma 2 q t - a 2 qTq ,(0) - a q T q (0)) 

where <j>\ is the master pre-stabilized laser phase. 

In the special case where the total round trip delay- 
times are equal, [T2 +T2>] = [t 3 +T3'], and the differential 
SC velocities are zero, P2 — — 0, the difference of the 
sensor signals, 



X 



(6) 



generates a standard equal-arm Michelson interferometer 
output, independent of laser phase noise. |23| For LISA, 
this is rarely a reasonable laser phase cancellation tech- 
nique since the arm-lengths are almost always un-equal. 
However, the TDI-Ai combination f24j, written as 



X x = A 2 - A 3 - A 2; 



22' 



*3;33' 



(7) 



replicates an equal-arm interferometer phase delay and 
cancels the common laser phase noise in the case where 
[T2 +T2>] ^ [t 3 + t 3 /] and /3 2 - /? 3 — 0. Calculating the 
timing-error between the </>i ;33 '22' and c/>i ; 22'33' terms in 
Eq.[7j which result from laser phase transformation order 
of the A 2; 22' and A 3;33 / terms, Eq.[TJ we find: 



8t = irlfa - p 3 \ 



(8) 



where r is the mean one-way delay time (~ 16.7s). Ex- 
ploiting Eq.QT] we can calculate the suppression limit of 
the TDI-Xi combination, which fails to account for this 
SC-motion induced timing error, as 22|: 



X x >4r|/3 2 -/3 3 |0i 



(9) 



where (f>i is the time-differentiated laser phase spectrum. 
Given a situation where this TDI-Ai limit is large enough 
to restrain the IMS sensitivity, the TDI-A2 combination, 



X2 =A2 — A 3 — A2 ; 22' + A 3;33 / 
— A2 ;33 '22' + A 3; 22'33' 

+ A2;22'22'33' — A 3;33 / 33 /22' 



(10) 



must be used to cancel the velocity coupled laser phase 
noise [20]. This TDI-X2 combination accounts for the 
timing error of the TDI-A"i combination by re-tracing 
and applying the TDI-Xi laser phase-delay transfer chain 
through the constellation a second time. 

The TDI-Xi and TDI-X2 combinations include multi- 
ple single-link, s sr , signals. Thus, the allowed noise for 
these TDI-Xi and TDI-X2 increase by a factor of 2 and 
4 respectively as compared to those referenced and ac- 
counted for in Table. [T] 

Although the following analysis will focus on the TDI- 
X2 velocity corrections, annual changes in j3, or SC ac- 
celeration terms, may be accounted for with further ex- 
pansion of these TDI combinations. Otherwise and in 
order to utilize the TDI-ranging methods outlined be- 
low for LISA TDI data-analysis, the f3 value will have 
to be adjusted to avoid the acceleration-dependent ac- 
cumulated error. Even though a continuous measure 
and correction to the f3 values is possible, these ad- 
justments may also be accomplished by segmenting the 
data-analysis and adjusting the ranging functions used to 
evaluate the TDI combinations at regular intervals, in a 
worst case LISA-like scenario, every 5t T year /(j3 it) — 
^/3.3 ns * 3.15 x 10 7 s/(66 ns/s * tt) = 708s [|. 



TDI Ranging 

Until now, it has been assumed that we know the re- 
quired interferometer arm-lengths and rate of change in 
order to form the TDI combinations but in practice this 
is not the case. Two primary methods, pseudo-random 
noise (PEN) ranging [H, and TDI-ranging 
have been proposed to measure these time-dependent 
arm-lengths. Extending the root mean squared (EMS) 
power minimization TDI-ranging methods outlined by 
Tinto in [14j . this experiment will, instead, modulate the 
laser phase signals with ranging reference tones at fre- 
quencies outside of the LISA measurement band. The 
EMS power minimization around these relatively high- 
frequency tones avoids the displacement of the measured 
arm-lengths as a result of low-frequency GW signals [l4[ 
and provides an improved ranging precision beyond the 
inherent EMS laser noise cancellation resulting from the 
larger signal power at these chosen frequencies. 

Using the Taylor approximation, 



X E rr{t) (f>(t)-<t>(t+5T), 



XErr{u) 
\X E rr(u)\ 



D — iujt _ — 2u;(i+5r)i 



u6t\4>\, 



(ii) 



for uiSt << 1, we can estimate a simplified but reason- 
able measure of the arm-lengths through the cancellation 
of these ranging tones using the TDI combinations to an 
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TABLE II. TDI generations based on orbital dynamics approximations || [2(j 



Generation Michelson Arm-Length Counter-Propagating Delay Delay Dynamics 



TDI-Ao.o 


T~22 


(*) = r 3 > 3 (t) 


T,(0) = 


v(o) 


dT q (t)/dt 


= 


TDI-Xi.o (First Generation TDI) 


T'22 


(*) + T 3 , 3 (t) 


T,(0) = 


v(o) 


dT q (t)/dt 


= 


TDI-A1.5 (Modified TDI) 


T 22 


(t) + T 3 , 3 (t) 


tM + 


v(o) 


drjt)/dt 


= 


TDI-A2.0 (Second Generation TDI) 


T22 


(t) + r s , 3 (t) 


tM * 


v(o) 


dr q (t)/dt -- 


= & 


TDI-X3.04. 


T22 


(t)^r 3 , 3 (t) 


T,(0) / 


v(°) 


d Tq (t)/dt = 


&(*) 




FIG. 2. Flow chart of the ranging-tone minimization process: This process minimizes the ranging tone and maximally constrains 
the six variable light travel time delays resulting in an optimized TDI-A2 strain sensitivity. The results of this process for the 
different experimental configurations are in Tablc lTVl 



accuracy of: 



SL = Stc : 



2tt/ 7 



-G 



Tone 1 



(12) 



where frone is the ranging-tone modulation frequency 
(1 Hz for these experiments) and Grone is the tone sup- 
pression magnitude when evaluated with the TDI combi- 
nation. Generally, the cancellation of the local laser sig- 
nal's ranging tones, <p\{t), from the far s sr signals in the 
TDI-X combinations (Eq.[7l [T0|) constrains the one-way 
outgoing delay times, T%(t) and Ts(t), while the cancella- 
tion of the far laser signal's tones, (f>2(t) and (fo(t), from 
the local s sr signals constrains the incoming delay times, 
Ty(t) and T2'(t), respectively. Exploiting the phase- 
locking methods outlined in 27| such that si 2 — S13 ~ 
or, equivalently, 4>2 = 0i ; 3' and 4>$ = </>i ;2 ', we can trans- 
fer the local laser phase data and phase stability to the 
far lasers and constrain the round trip delay times: 



alt 



T 2 2'(t) — ul 2 
733' (t) = ajt 



o| T2 ,( ) 



a 2 T2(0), 
a2T 3 (0), 



(13) 
(14) 



using a local ranging tone only. 

Although some prior estimate of the arm-lengths will 
likely exist, the analysis in this paper will assume no pre- 
vious knowledge of the 6-variable time-dependent round- 
trip arm-lengths, T22'(t) and T33/(f), and will determine 
these arm length functions using a 1 Hz laser frequency 
modulation with an amplitude of 500 Hz. Note that the 
ranging tone should not be placed at a frequency which 
is near an integer multiple of the interferometer arms' 
inverse round-trip delay time to avoid the inherent tone 
cancellation along a single arm as shown in Fig.[3l This 
local ranging tone will provide constraints on the round 
trip delay times only. Since the far laser phase signals 



in this experiment are not modulated with ranging tones 
but rather, phase locked [27j], the one-way delay times 
are constrained using the TDI-ranging methods outlined 
by Tinto (Hf through the minimization of the phase lock 
loops's (PLL) residual phase RMS power in the TDI-X2 
combination. 

As shown in Fig. [21 initially assuming the arm-lengths 
are constant (/3 = 0) a four-dimensional sweep of the 
'time-space' defined by r 2 (0), r 2 ' (0), T3(0), and r 3 '(0) is 
performed using a 51-point Lagrange fractional delay fil- 
ter [28} which determines the values which minimize RMS 
power near the 1 Hz frequency-modulated ranging tone 
in the TDI-X 2 combination. Generally, depending on 
the pre-stabilized laser noise and phase locking configu- 
ration, the delay times must be measured with a 10 -8 
resolution 0, [23. Therefore, to make data analysis more 
efficient, rather than evaluating the 10 s possible values 
along each of the 4 delay-dimensions, the time-delay de- 
termination is performed using an algorithmic scan with 
successively finer time-delay mesh grids until the ranging 
tone is minimized and dominated by the noise floor of 
the IMS measurement. 

Using this scanning process, we determine an ini- 
tial measure of the round-trip delay offsets, r 22 /(0) and 
T33'(0), which minimize the ranging tone for small data- 
segments along a continuous data-set. The fitted slope 
of these offsets evaluates a first-order approximation of 
/3 and the time-dependent arm-lengths. Applying the j3- 
value dependent time-expansion or time-contraction to 
the s sr signals using time- varying fractional delay inter- 
polation [28[, the process is repeated iteratively, further 
optimizing the arm-length functions as shown in Fig. [2] 
Finally, the fitted values, /3 2 , f3 3l t 2 (0), t 2 / (0), T3(0), and 
T3'(0), are used to evaluate the TDI-A 2 combination 
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FIG. 3. Spectral magnitude tone suppression: The ranging tone modulated inputs are plotted along with the resulting TDI-A^ 
combinations for the four different experiments outlined in Table lTTTl The decreased round-trip timing accuracy of the A2 arm 
as compared to the A3 arm as shown in Table lfVl is likely due to the proximity of the nearest arm-zeros. 



along the entire data-set. The results produce a measure 
of the round-trip arm-lengths and place constraints on 
the one-way arm-lengths to an accuracy beyond the rang- 
ing requirements, thus, removing the sensitivity-limiting 
laser and PLL phase noises sources and idealizing the 
total interferometer strain precision. 



LISA SIMULATOR BENCHTOP 

The University of Florida LISA Simulator (UFLIS) 
benchtop (Fig.0} consists of four lasers, three of which, 
Laser-1 (Li), Laser-2 (L2), and Laser-3 (L3) represent the 
lasers on each of the SC in the LISA constellation. Laser 
beatnotes are formed between each of these laser fields 
and a global reference laser (L^), acting as an optical 
clock against which each individual laser phase is mea- 
sured. Combinations of these PD measurements cancel 
the common Lr phase noise. Lr and Li are stabilized 
through Pound-Drever-Hall (PDH) stabilization [29| us- 
ing ultra-low expansion reference cavities to achieve a 
~ 100 Hz/vTIz pre-stabilized laser noise as shown in 
Fig. [6] Even though a lower laser frequency noise is 
achievable it is intentionally spoiled to display the 
suppression capabilities of the TDI combinations. 



Phasemeter 

The phasemeter is used to measure the phase of 
a 2-20 MHz PD beatnote signal to an accuracy of ~ 
1 /j,cycle/-\/Hz (Table|T|. The beatnote is sampled using 
a 14-bit analog-to-digital Converter (ADC) at a rate of 
40 MHz. The sampled signal is processed by a field pro- 
grammable gate array (FPGA) programmed with a dig- 
ital PLL. The digital PLL's feedback signal is down- 



sampled to a rate of 19.1 Hz and relayed to a data- 
processing computer. Using differential and entangled 
phase |l3L l30j PM measurements it has been determined 
that the PM is limited by <Padc (FiglU, a combina- 
tion of frequency-dependent ADC timing jitter (5t = 
35 / V7 fs / VHz) j RF-transformer phase dispersion, and 
amplitude noise [30j . One may write the PM measure- 
ment output as 

fin 

4>PM{t) =<f>in(t) + T^— 4>Clk(t) + 4>ADc{t), (15) 
JClk 

where <f>i n is the phase information on the /i„-frequency 
beat-note, 4>cik is the phase noise of the /czfc-frequcncy 
sampling clock, and 4>adc are the ADC-noise sources 
mentioned above. Note the coupling of the clock's phase 
noise into the measurement. PM measurements taken on 
different LISA SC will be clocked using different ultra- 
stable clock sources requiring the need for clock-noise 
transfers between SC to remove these terms [31]. Al- 
though the following TDI experiments are taken use a 
single clock, the TDI combinations will still require clock- 
noise corrections to account for the time-shifted PM clock 
noise terms. This is discussed in the following section. 

Electronic Phase Delay Unit 

The EPD unit simulates the characteristics of the laser 
phase transmission between the SC including the time- 
varying light travel time, Doppler shift, and gravitational 
wave phase modulations using a high-speed digital signal- 
processing (DSP) system. The front end is a fast PM, 
similar to the PM described above, with a data rate 
of 19.53 kHz instead of 19.1 Hz. The PM data is in- 
terpolated to time-lead or time-lag the phase informa- 
tion, producing a linear variation in the time-delay. GW 
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FIG. 4. Experimental Model of the LISA Interferometry Benchtop: The arrangement of the laser phase sources, EPD units, 
phase-lock loops, PD signal beatnotes, and PM measurements which are used to emulate the LISA-IMS in our experiments are 
presented. The measured s sr signals are used to form the TDI combinations. 

TABLE III. TDI Experimental Characteristics: Four 40000 s experiments are performed with increasingly more complicated, 
yet more LISA-like, characteristics. The static transponder experiment provides us with a baseline measure of the experimental 
setup's noise performance. The dynamic transponder experiment demonstrates the ability to determine and account for the 
time-changing delay-times. The phase-locked LISA-like experiment proves the ability to remove independent SC noise sources 
and constrain one-way delay times. Finally, the confusion noise experiment verifies that the TDI-ranging capability will not be 
limited by low-frequency LISA noise sources. 
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signals and a Doppler offset are added to interpolated 
phase information before it is used to drive a numer- 
ically controlled oscillator (NCO). Completing the laser 
phase transmission replication, the NCO output is regen- 
erated using a 16-bit digital-to-analog converter (DAC) 
with the same clock-source as the PM. After accounting 
for the clock noise coupling and Doppler shifts, we can 
write the EPD unit's output as, 

4>BPD{t)=4> in {t-T{t)) (16) 
+ <t>Anc{t - r(t)) + (f) DAC (t) 
+ 4>Clk:ADc{t - r(i)) - 4>Clk:DAc(t) 

+ fin f fDOP [<j>Cl k (t-T(t))-<f>c lk (t)}- 
JClk 

4>in(t) is the phase evolution of the beatnote at an av- 
erage frequency of /j„ = 2 — 20 MHz. The beat signals 



are measured and regenerated with respect to a single 
clock, 4>cik(t), at the DSP system clock frequency of 
fcik — 40 MHz. The phase noise of this clock enters as 
phase variations between the sampling and the Doppler 
shifted, /_Do P , regeneration after a time- varying time- 
delay r(t) = r(0)+j3t. The single clock source is split and 
distributed between the sampling ADC and regeneration 
DAC causing an additional phase error, (f>cik-.ADc(t — 
f{t)) - (f>cik-.DAc(t). The sampling (4>adc) and regen- 
eration (4>dac) processes add additional converter-noise 
contributions. It has been determined using two different 
sampling and regeneration clock sources that the limit- 
ing noise source on the EPD unit (FigJSJ) is the differen- 
tial clock phase noise, (j>cik:ADC,r(t) - <t>cik-.DAC, caused 
by errors in the clock distribution [301 ] - Note, as one 
may check, that the clock noise terms, cj>cik{t) in Eq.fTHl 
themselves will cancel when measured with phasemeters 
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fEq.[T5j) and evaluated in the TDI combinations (Eq.[7l 

ma. 



Experimental Setup 

The laser benchtop, PMs, and EPD units are combined 
to recreate the LISA-IMS (Fig.!}. The Lx/Lr differen- 
tial beatnote phase represents the pre-stabilized 'input' 
noise. This PD signal is electronically mixed with a 1 Hz 
frequency-modulated oscillator to add a ranging-tone and 
produce the input laser phase signals. Replicating an 
interferometer beam-splitter, these signals are electroni- 
cally split and processed by the EPD units to simulate 
the outgoing-held inter-SC light transmission. The EPD- 
processed signal is mixed with each of the L 2 / 3 /Lfl beat- 
notes 'on the far spacecraft' to produce the S12 and S13 
PM signals. These delayed signals are also used to phase- 
lock L2/3, transferring the Li stability to these lasers 
[27] ]. The L 2 / 3 /Li? beatnotes are again processed by EPD 
units, simulating the returned-field inter-SC light trans- 
mission. Finally, the local differential Li/L/j input phase 
signal is mixed with the delayed L2/3/LA beatnotes to 
produce the s 2 i and s 3 i PM signals. It should be ex- 
plicitly noted that these LISA-like s sr signals only in- 
clude noise due to the pre-stabilized laser noise source, 
the phase-lock loop's phase noise, and EPD unit's abil- 
ity to mimic the inter-SC time delay and are not sensi- 
tive to path-length differences or the shot-noise limit of 
the laser bench itself. Excluding these displacement and 
shot noises from the IMS noise budget reduces the single 
link requirement to about 5pm/vTIz, the TDI-A"i com- 
bination requirement to 14.1 pm/\/Hz, and the TDI-A 2 
combination requirement to 20 pm/VHz 

Implementing this experimental model, four measure- 
ments are performed as outlined in Table lHTl In the 
transponder measurements, rather than phase-locking 
the outgoing-field's EPD signal, it is transferred directly 
to the return-field's EPD unit as though it were re- 
flected off of a moving mirror; accordingly, these TDI 
experiments may completely neglect the s± r terms while 
ranging is only required for the round-trip delay times, 
r 22 '(t) (Eq.HS]) and r 33 <(t) (Eq.E]), instead of all four 
one-way delay times. For all measurements presented, 
the arm-lengths are chosen as r 2 ~ r 2 ' ~ 16.55 s and 
T3 ~ T3< ~ 16.75 s to maximize the unequal arm-length 
mismatch. The relative spacecraft velocities for the dif- 
ferent measurements are shown in Table lHTl They are 
artificially large to increase the |/3 2 — /?3 1 limitations (Eq. 
Ej]) for the TDI-Xi combination and to prove the ability of 
the TDI-A 2 to remove the differential velocity dependent 
noise from the TDI-Ai combination. Doppler shifts of - 
2.0 MHz and +3.0 MHz are applied in order to produce 
the necessary MHz beatnote observables. A frequency 
modulated 6.22 mHz verification binary GW signal with 
an amplitude of 5 few — 2 /iHz is injected into all four 




FIG. 5. Static Transponder (Baseline) Experimental Re- 
sults: The sensor signals (s2i(i) = ^{t), S3i(t) = As(t)) 
are plotted along with the raw TDI-Xi and TDI-X2 results. 
The phasemeter measurement limitation and expected rang- 
ing limitations based on the calculated timing variance are 
also plotted. 




FIG. 6. Corrected Static Transponder (Baseline) Experimen- 
tal Results: The input signal, <f>\{t), is plotted together with 
the TDI-Xi and TDI-X2 results, corrected by their respective 
transfer functions. Both TDI combinations are limited by the 
noise added by the EPD unit (grey curve). 



measurements to verify GW extraction. This frequency 
modulation equates to a phase modulation amplitude of 
S^GW = <5 few I (27T '/) = 51.2/zcycles resulting in a one- 
way strain amplitude of 5h — 5(/>[X/(ct)] = 1.1 x 10 -20 . 
The resulting GW strain amplitude h = A8h = 4.4 x 
10~ 20 is a factor of 100 larger than expected for the RX 
J0806.3+1527 AM CVn binary [32]. Finally, a mock- 
confusion noise is added to the signals to demonstrate 
that this low-frequency noise will not limit the ranging 
capabilities. 

The following experimental results are averaged over 
six, 10000 s, data-segments during the course of a contin- 
uous 40000 s experimental run-time. 
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FIG. 7. Laser Noise Magnitude Suppression Function: The 
achieved transponder-measurement laser noise suppression 
magnitudes of the TDI-Xi and TDI-X2 combinations are 
plotted for both the static and dynamic experiments. The 
TDI-Xi combination's noise suppression equals the theorized 
limit (Eq.EJ 



TDI RESULTS 

Static Arm Transponder 
(Baseline) 

We first re-establish [13j a baseline ranging and mea- 
surement precision with static arm-lengths (/3 — 0), uti- 
lizing the PD measurements of S21 and S31 and assum- 
ing Si2,si3 — when evaluating the TDI combinations. 
The 40000 s data-set is broken into 40, 1000 s segments. 
The first iteration and linear regression of the ranging 
process described in Fig. [2] produces a slope error (con- 
straint on the arm-length velocities), of \2f3\ < 50fs/s 
and a variance (round-trip ranging accuracy) of 0.6 ns 
(~0.18m) as shown in Table llVl In this experiment, we 
note that the TDI-Xi combination's ranging-tone mini- 
mization produces the same result as the TDI-X2 com- 
bination to within the measurement error. If /3 ^ 0, this 
would not be the case since the ranging tone minimiza- 
tion using the TDI-Xi combination would be limited by 
Eq.[5] and would tend to calculate the mean delay for a 
particular data- segment. 

Using the calculated values we form the TDl-Xi and 
TDI-X2 combinations for the entire data-set. The raw 
TDI results, as plotted in Fig.[5l show the laser noise 
cancellation and reveal the phase-modulated GW binary 
at 6.22 mrlz. The residual noise is dominated by phase 
noise added by the EPD units and, based on the timing- 
error, not by uncertainties in the ranging. The differ- 
ences between the TDI-Xi and TDI-A2 combinations' 
spectral noise are caused by differences in their trans- 
fer functions with respect to the input laser phase noise. 
Fig.[6]shows the spectra after the TDI combinations have 
been rescaled by their respective transfer function magni- 
tude. Both agree well with each-other and demonstrate 



greater than 10 orders of magnitude laser phase noise 
suppression below 1 mHz. (Fig. [7]) 

Dynamic Arm Transponder 
(TDI 2.0 Verification) 

In the next experiment, we incorporate time- 
dependent arm-lengths into the simulation with the (3- 
values defined in Table lTTTl Again, initially assuming f3 = 
0, the 40000 s measurements of the S21 and s 3 i signals are 
broken into 40, 1000 s segments. These data segments 
are then used to minimize the ranging tone and calculate 
the round-trip time-delay for each segment as defined by 
Fig. [2] The linear regression of these time-delay offsets 
provides a first-order measure of the (3 to an accuracy of 
lOOps/s as shown in TablellVl The process also calculates 
a first order measure of the round-trip time delay with a 
ranging precision of < 7.5 /is (^1.7km) but, due to the 
incorrect /3 = assumption, these values tend to equal 
the average delay for the data-segment. A second itera- 
tion improves the (3 accuracy to 80fs/s and the ranging 
precision to < 5.9 ns (~1.8m). The final iteration opti- 
mizes the (3 precision to 8.9 fs/s and the ranging precision 
to < 0.65 ns (~0.2m). (Table lTV| Producing comparable 
values to within the measurement error, additional itera- 
tions do little to improve the tone cancellation or ranging 
accuracy. 

Applying the calculated round-trip functional values 
from the third iteration of the ranging procedure, we use 
the entire data-set to produce the TDI-JTi and TDI-X2 
combinations (Fig. [8]). The TDI-JTi combination is lim- 
ited, as theoretically anticipated, by Eq.|H]with t ~ 16.7 s 
and I A/3 1 = 250ns/s. The TDI-X2 combination's cor- 
rection terms account for this dynamic arm-length lim- 
itation and remove the velocity dependent laser phase 
noise resulting in a sensitivity equal to the experiment's 
baseline noise performance. This result meets the IMS 
requirement defined by the LISA mission concept design. 
(Table[T| The ranging precision, as plotted in Fig.[8l is 
not expected to be a limiting noise source which is veri- 
fied with through the cross-correlation of the input noise 
with the TDI-X2 combination as shown in Fig.HUl 



Dynamic LISA-like 
(Phase-locked Laser, One-Way Ranging) 

At this point, we include the phase- locking of the L 2 /3 
signals on the far spacecraft (Fig.HJ) and the transmis- 
sion of these phase signals back to the local-SC, thus 
generating and measuring all four s sr beatnote observ- 
ables. These signals are used in the same iterative pro- 
cess as previously described (Fig. [21 Table lIV|) . The op- 
timized time-delay functions from this process result in 
a measure of \2f3\ to an accuracy better than ^70 fs/s 
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TABLE IV. Ranging Precision 



TDI Experiment Name 












TDI-Ranging Constraint 












Iteration (TDI combination) 


a 
P 




T 22 i (V ) , T33' (U ) 


OT 2 2> 


> T 33' 


Static Transponder 












Round-trip Ranging 












1 (TDI 1.0) 


2/3 2 = -44.5 fs/s 
2f3 3 = -46.3 fs/s 


±20.9 fs/s 
±12.5 fs/s 


r 22 ,(0) = 33.55204887148 s ±0.23 ns 
733' (0) = 33.15222859583 s ±0.14 ns 


St 22 , 
St 33 > 


= 0.54 ns 
= 0.32 ns 


1 (TDI 2.0) 


2/3 2 = -41.0 fs/s 
2/3 3 = -46.3 fs/s 


±21.2 fs/s 
±12.7fs/s 


t 22 /(0) = 33.55204887151 s ±0.24 ns 
733,(0) = 33.15222859579 s ±0.14ns 


5t 22 , 
St 33 > 


= 0.55 ns 
= 0.33 ns 



Dynamic Transponder 



Round-trip Ranging 



1 (TDI 2.0) 


2/3 2 
2f3 3 


= -200.247 ns/s ±100ps/s 
= ±300.056 ns/s ±95ps/s 


T22'(0) 
7-33' (0) 


= 33.5518847 s ±2.3 ^s 
= 33.1525027 s ±2.2 ps 


8t 22 > 
5t 33 > 


— 7.5 fis 
= 7.0 fis 


2 (TDI 2.0) 


2f3 2 
2f3 3 


= -199.9998668 ns/s ±80 fs/s 
= ±300.0001130 ns/s ±23 fs/s 


722' (0) 
T 33 >(0) 


= 33.5519484187 s ±1.8 ns 
= 33.1523897572 s ±0.51 ns 


8t 22 > 
5t 33 > 


= 5.9 ns 
= 1.7 ns 


3 (TDI 2.0) 


2f3 2 
2/3 3 


= -200.0000058 ns/s ±8.9 fs/s 
= ±300.0001361 ns/s ±4.5 fs/s 


T22<(0) 
7-33' (0) 


= 33.55194832884 s ±0.20 ns 
= 33.15238977691 s ±0.10 ns 


8t 22 , 
5t 33 i 


= 0.65 ns 
= 0.33 ns 


Dynamic LISA-like 














Round-trip Ranging 














1 (TDI 2.0) 


2/3 2 
2/3 3 


= -199.984 ns/s ±12ps/s 
= ±300.052 ns/s ±7.8ps/s 


T 22 '(0) 
7-33' (0) 


= 33.59821021s ±0.28 ^s 
= 33.21476669 s ±0.18^s 


5t 22 , 
St 33 > 


= 0.895 ^s 
= 0.568 ^s 


2 (TDI 2.0) 


2f3 2 
2/3 3 


= -200.000015 ns/s ±71 fs/s 
= ±300.000013 ns/s ±26 fs/s 


7-22' (0) 
7-33' (0) 


= 33.5982645303 s ±1.6 ns 
= 33.2146434958 s ±0.58 ns 


5t 22 , 
St 33 > 


= 5.2 ns 
= 1.9 ns 


3 (TDI 2.0) 


2P 2 
2/3 3 


= -200.000028 ns/s ±68 fs/s 
= ±300.000020 ns/s ±25 fs/s 


7-22' (0) 
7-33' (0) 


= 33.5982645401s ±1.5 ns 
= 33.2146435166 s ±0.58 ns 


8t 22 > 
St 33 > 


= 5.0 ns 
= 1.9 ns 


One- Way Ranging 


3 (TDI 2.0) 


02 

03 


= -103.3 ns/s ±1.4 ns/s 
= ±152.24 ns/s ±1.4 ns/s 


7-2(0) 
7- 2 * (0) 
7-3(0) 
7-3*' (0) 


= 16.68021s ±32 ^s 
= 16.91805 s ±32 fjs 
= 16.48824 s ±31 fis 
= 16.72640 s ±31 ^s 


St 2 
5t 2 > 

5T3 

St 3 > 


= 105 fis 
= 105 fis 
= 99 fis 
= 99 fis 


Dynamic LISA-like 
with Confusion Noise 














Round-trip Ranging 














3 (TDI 2.0) 


2f3 2 
2/3 3 


- -199.999991 ns/s ±50fs/s 
= ±300.000137 ns/s ±22 fs/s 


7-22' (0) 
7-33' (0) 


= 33.6012734891 s ±1.1 ns 
= 33.2100302983 s ±0.49 ns 


5t 22 , 
Sr 33 > 


= 3.7 ns 
= 1.6 ns 


One- Way Ranging 


3 (TDI 2.0) t 


02 
03 


= -96.81 ns/s ±2.3 ns/s 
= ±149.439 ns/s ±1.4 ns/s 


7-2(0) 
7-2* (0) 
7-3(0) 

7-3* (0) 


= 16.73546 s ±53 fis 
= 16.86582 s ±53 fis 
= 16.53994 s ±32 fis 
= 16.67009 s ±32 fis 


St 2 
5t 2 > 
5t 3 
St 3 , 


= 169 fis 
= 169 fis 
= 102 fis 
= 102 fis 



J The additional error in the one-way confusion noise measurement as compared to the phase-locked measurement is 
expected to be due to the coupling of low- frequency noise using the minimized-RMS ranging method [l3 |. 



* The returning delay times tend to be longer than the outgoing delay times by 100-300 ms as a result of internal delays 
within the DSP system's EPD units. 
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FIG. 8. Dynamic Transponder Experimental Results: The 
suppression of the TDI Xi combination is limited by the arm- 
length time-dependence. The TDI X% combination removes 
the additional time-dependent-coupled laser noise and reveals 
the 6.22 mHz GW signal. As with the static case (Fig. [6]), 
the EPD unit's phase transmission accuracy is the primary 
limiting noise source in the TDI combinations although, some 
sensitivity loss may occur due to a limited ranging capability 
around 100 mHz (Fig.H0|. 



and a round-trip ranging precision of ^5. 0ns (~1.5m). 
The constraints on the one-way delay times through the 
residual PLL noise removal (~1 mHz/VHz) are not ap- 
plicable until the precision of the round-trip delay times 
are accurate enough to remove enough of the input laser 
noise from the TDI combinations to reveal these residual 
PLL noises. Thus, it requires at least one iteration of the 
ranging process before one can constrain the one-way de- 
lay times. Applying a linear regression to the calculated 
one-way delay times we find a one-way ranging error of 
~ 100 /is (^30 km). The outgoing and return delay times 
are un-equal by ~ 250 ±0.1 ms, proving the ability to ex- 
tract the individual one-way laser phase errors despite 
un-equal delays along a single arm (r 9 (0) ^ T q >(0)). 

Applying these optimized one-way functional values 
from the ranging procedure, we produce the TDI-Xi 
and TDI-A" 2 combinations (Fig.|9j. Again, the TDI-Xl 
combination equals the expected limitation (Eq.|Hl). The 
TDI-X2 combination meets the LISA IMS requirement to 
within a factor of 4 and is likely limited by a combination 
of multiple EPD clock-noise sources resulting in a sensi- 
tivity greater than the simulator's baseline performance. 
Based on the variance of the fitted delay times, the rang- 
ing precisions are not a limiting noise source as plotted in 
Fig-El The cross-correlation of the TDI-X2 combination 
with the laser and PLL noise sources (Fig.[T0j) indicates 
that all the known and accountcd-for noise sources have 
been sufficiently removed from the interferometer's out- 
put. 



Dynamic LISA-like 



FIG. 9. Dynamic LISA-like (Phase-locked) Experimental Re- 
sults: The suppression of the TDI Xi combination is limited 
by the expected arm-length time-dependence. The TDI X2 
combination removes the input laser noise, the far-SC PLL 
residual phase noise, and the time-dependent coupled laser 
noise to reveal the 6.22 mHz GW signal. The sensitivity limi- 
tation comes, most likely, as a result of multiple uncorrelated 
EPD noise sources. 
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FIG. 10. TDI vs. Input Cross Correlation: The magni- 
tude squared cross-correlation of the LISA-like TDI measure- 
ments show no correlation with the input noise or either PLL 
noise source, verifying that all the laser noise sources have 
been sufficiently removed. The cross-correlation of the dy- 
namic transponder TDI measurement shows some input phase 
correlated noise cancellation limitation for frequencies above 
100 mHz. 



(with Confusion Noise Background) 

Lastly, we include a low frequency simulated 'con- 
fusion noise' into the measurement to ensure that 
these low-frequency terms do not limit the rang- 
ing precision. The confusion noise background, 
6.4/(A/l + U /Ir) 2 ) /icycles/VSi where 
32 



'{fy/l + (/// fl ) 2 )/xcycles/v / Hz where f R = 1 mHz 
, and the 6.22 mHz mono-chromatic binary are simul- 
taneously injected. The optimized ranging result places 
bounds on the \2/3\ accuracy better than ^50fs/s pro- 
duces a round-trip ranging precision of ~3. 7ns (~1.1 m). 
Thus, this confusion noise result achieves a ranging pre- 
cision on the same order as the simulator's phase-locked 
performance, indicating that low-frequency noise has lit- 
tle to no effect on the ranging tone cancellation or the 
measured arm-lengths. 

The noise spectrum comparisons between the TDI- 
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X2 outputs of the dynamic LISA-like experiments, with 
and without the confusion noise background, are plot- 
ted in Fig.QTJto show the additional low- frequency noise. 
These spectra are a factor of 5 larger than the simula- 
tors baseline performance likely resulting from the cou- 
pling of the additional EPD, phasemeter, and electronic 
components in these experiments. All three measured 
spectra in Fig. [11] have been scaled by the high-frequency 
sensitivity loss 'roll-up' of the LISA-detector for GW- 
frequencies larger than 1/r = 60 mHz in order to obtain a 
direct comparison with the single-link sensitivity. The in- 
jected confusion noise background is plotted and matches 
the spectrum of measured confusion noise. A theoretical 
confusion-noise background |33j and the expected 1-year 
strain amplitude of the four strongest LISA verification 
binaries 32| are marked for additional reference. 

The time-series of the extracted monochromatic GW 
binary signal using the TDI-A2 combination is plotted in 
Fig.[l2]with and without the confusion noise background. 
Comparing the measured amplitude with the expected 
GW amplitude, given the injected 51.2 //cycle EPD GW- 
signal, we find the TDI-A2 extracted signal matches the 
expected amplitude of 4 x 51.2 /icycles= 205 ^cycles. The 
factor of 2 accounts for the two interferometer arms while 
the factor of 2.33 accounts for the TDI transfer function's 
signal gain at / = 6.22 mHz (|A 2 (6.22 mHz)| = 2.33). 



CONCLUSION 

Following the initial interferometry tests of a static 
LISA model [Hj], we expanded UFLIS and added time 
varying signal travel times, Doppler shifts, and gravita- 
tional wave signals to our electronic phase delay units. 
This enables tests of the LISA interferometry in a realis- 
tic, dynamic model. Our experimental results show that 
than 10 orders of magnitude of laser phase noise can be 
canceled using appropriately time-shifted data streams in 
the TDI-A 2 data combination. We also developed and 
demonstrated a simple but powerful ranging method to 
measure the signal travel times between the spacecraft. 

We verified that the ability to reduce laser phase noise 
using a TDI-Ai data combination is indeed limited by 
the relative velocities between the spacecraft. Further- 
more, we demonstrated the removal of the residual phase 
lock loop noise added at the far spacecraft and, in this 
configuration, showed that the requirements on one-way 
ranging are relaxed by several orders of magnitude com- 
pared to the requirements on round trip ranging. 

In the process, we developed and tested data analysis 
tools which take the raw phasemeter data streams, ex- 
tracts the light-travel time functions, and generates the 
TDI-A 2 data sets. We also added a confusion noise GW- 
background to our interferometry emulator and verified 
that this background noise does not interfere with our 
ranging capabilities. 
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FIG. 11. Compiled Results and Comparison with TDI for 
LISA: In this figure we have compiled all the results of the 
TDI simulations and attempt to make a direct comparison 
with the expected LISA strain sensitivity. The baseline spec- 
tral noise of the UFLIS simulator (grey-blue) from the TDI- 
Transponder measurements is plotted. The velocity corrected 
TDI-A2.0 spectrum of the dynamic arm TDI simulation with 
(cyan) and without (blue) the injected binary confusion noise 
(dotted-magenta) is plotted in comparison with the IMS sen- 
sitivity requirement. The three TDI simulations are scaled 
to account for the high-frequency GW-sensitivity loss ex- 
pected in LISA. The DRS and IMS requirement are root- 
square summed and scaled by the high-frequency LISA GW- 
sensitivity loss function to produce the effective single-link 
LISA sensitivity. An estimate of the confusion noise limit is 
plotted (dotted-red) along with the four brightest verification 
binaries rescaled from a 1-year averaged strain sensitivity to 
noise spectra in cycles /V Hz. The strain magnitude of the 
1 year averaged RX-J0806 binary and the 10000 s EPD in- 
jected GW have amplitudes such that they result in similar 
LSD amplitudes in this figure. 



Measured TDI-2 withMonochromatic Binary + Confusion Noisi 
- Measured TDI-2 with Monochromatic Binary 
■ Expected TDI-2 GW Amplitude 




FIG. 12. TDI-X2 Time-Series: The TDI-X2 extracted signals 
match the expected GW amplitude given the 51.2 /icycle GW 
EPD injection amplitude after scaling by the magnitude of 
the TDI-X2 combination's transfer function evaluated at / = 
6.22 mHz, \X 2 (6.22 mHz) | = 2.33. 



Future experiments should include real, LISA-like GW 
signals using data-sets generated with LISA-tools like 
Synthetic LISA [sj. Simulations with three indepen- 
dently stabilized lasers might also be valuable towards 
verifying the constrains on the one-way ranging capabil- 
ities. 
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